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^^ I We study the role of finite widths of resonances in a nonlocal version of the Wick-Cutkosky 

f^ , model. The spectrum of bound states is known analytically in this model and forms linear Regge 

f^ ' tragectories. We compute the widths of resonances, calculate the spectral density in an extension of 

CN| ' the Breit-Wigner ansatz and discuss a mechanism for the damping of unphysical exponential growth 

.^.j , of observables at high energy due to finite widths of resonances. 

o . 

Q ' PACS numbers: 12.38.Aw 12.38.Lg ll.15.Tk 11.55Bq 11.55 Jy 

^n ; 

'"' ; I. INTRODUCTION 

(N ; 

^ \ Behaviour in field theories with nonlocal interactions in the time-like momentum region is the central concern of 

0^ • this paper. Nonlocal quantum field theory has a long history. Initially it was invented as an attempt at solving the 
^^ ' problem of UV divergences. However after the development of renormalisation techniques for local QFT the study 
St . of nonlocal interactions became a field of academic interest. Nevertheless the challenge of constructing self-consistent 
^-v ■ and mathematically rigorous quantum field theories with nonlocal interactions and their phenomenological application 
■r-l- , have been attracting permanent attention. As a result of this long activity of several decades the general principles of 
f^ • quantisation of nonlocal field theories have been formulated. Various issues related to quantisation, unitarity of the 
S'-matrix, realisation of causality, the validity of Froissart-type bounds at hi gh e nergy and the physical interpretation 
of nonlocal fields can be found for instance in [llQiSllilSSHIaiS H^ Hll ■ For the present paper the essential 

i ' result is that for a theory with typical growth of amplitudes in the complex momentum plane and A the scale of 

^-1. nonlocality the Froissart type bound on the total cross-section •4'j , 
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atot{s)<Cs^-^\og^s, (1) 

^~e(p'/'^')^ K>l/2 
can be derived assuming 

• unitarity of the S'-matrix, 

• analyticity of amplitudes on the Martin-Lehmann ellipse. 

These two requirements are fulfilled in local as well as in nonlocal QFT with the form factors being entire analytical 
functions of momenta. Despite this upper bound on the total-cross section, practically useful procedures for calculating 
observables at high energy in nonlocal QFT have been not formulated. Thus, since amplitudes are exponentially 
growing within perturbation theory, the common wisdom is that nonlocal theories necessarily display patholo gica l 
behaviour at high energy, which is usually considered as a fundamental drawback of nonlocal QFT models |lll Il2| . 
The aim of the present paper is to give an example of such a practical calculation for the total annihilation cross 
section in a confining nonlocal model exhibiting a Regge spectrum of bound states. In order to decode this statement 
of the problem we first have to explain its context. 
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An essential disadvantage of nonlocal theories has always been seen in the functional ambiguity in the choice of form 
factors related to the nonlocal character of interactions. Though requirement of unitarity restricts the form factors 
to be entire analytical functions of momenta, an ambiguity within this class of form factors is unavoidable as far as 
nonlocality is introduced as a fundamental property of the theory. An unexpected source for a nonlocal character of 
quantum fields has been pointed out by Leutwyler 7\ who has noticed that the propagator of the massless charged 
scalar field in the presence of background (anti-)self-dual homogeneous gauge field takes the form 

D{p') = '—^, , (2) 

here p is Euclidean momentum, and scale A is related to the strength of the background gauge field. The nonlocal 
form factor (1 — exp(— p^/A^)) appears here due to the presence of strong nonperturbative gauge fields, and thus 
its appearance itself and particular form can be related to the properties of the physical vacuum of the system. For 
vanishing field strength A — > 0, and the usual scalar propagator is recovered. Such an observation was immediately 
recognised as potentially interesting in application to the confinement problem in QCD 7] . For physical interpretation 
it is important that the nonperturbative background gauge field eliminates the pole in the momentum space propagator 
of a charged field rendering it an entire analytical function. Such an analytical property of the propagator can be 
regarded as dynamical confinement of the charged field since it has no particle interpretation [3, 0, [l3 ■ It is also 
important that such well established perturbative properties of QCD as asymptotic freedom at short distances are 
unaffected by the nonlocality generated by such background gluon fields. We have here an example of a theory defined 
by a local classical Lagrangian but with nonlocal behaviour of the quantised fields appearing as a consequence of a 
nonperturbative gauge field configuration. The (anti-)self-dual homogeneous gluon field itself can be seen as a highly 
idealised analytically tractable realisation of more complicated and realistic nonperturbative gauge configurations 
responsible for confinement (for instance, the domain model [lj,ll5j is a specific attempt at realising a more detailed 
picture of the QCD vacuum in this spirit). 

Turning to confinement specifically now, this phenomenon is more than just the absence of the coloured quark- 
gluon states in the physical spectrum of QCD, but the mechanism by which colourless bound states are present in 
the spectrum of QCD. One of the first connections between confinement and the hadron spectrum was achieved by 
the string picture, which provided a simple explanation for the pattern of Regge trajectories from a mechanism of 
confinement. However, the string explanation for Regge trajectories is not unique. Another, purely field theoretical 
approach to bound state formation is provided by the hadronisation procedure in the QCD functional integral which 
is closely connected with Bethe-Salpeter approach to the description of relativistic bound states [la. Il7| . It has been 
shown within the bosonisation procedure that nonlocality in quark propagators similar to Eg . ^ generated by the 
homogeneous (anti-)self-dual field in QCD leads to asymptotically linear Regge trajectories 17] . In fact this approach 
turned out to be successful in application to wide variety of meson spectra - light, heavy-light mesons and heavy 
quarkonia |18J . The specific chiral properties of quark eigenmodes are similar to those in an instanton background 
and are of crucial importance in the description of light mesons in this context. 

In its most clean and refined form the role of nonlocality of the type Eq. |(2Jl in generating a linear Regge spectrum 
of bound states has been exposed in [13 where the method of [13 has been improved and applied to a scalar field 
model. Unlike genuine QCD with all its complications, this model allows for analytical solution of the Bethe-Salpeter 
equation in the one-boson exchange approximation. The Lagrangian of the model has the same structure as that of 
the well-known Wick-Cutkosky model (23 

£ = _$t5'-l(_a2)^ _ }_^j^-1^_q2^ _ g$t$0 (3) 

The Wick-Cutkosky model corresponds to a standard real massless scalar field and massive charged scalar field $. 
It was invented as a prototype theory for studying the relativistic bound state problem in quantum electrodynamics. 
It turned out that the choice of Euclidean momentum space propagators of the form 

S{p^)^D{p') = ^^e-^''^' (4) 

leads to a Bethe-Salpeter equation which is analytically soluble in the one-boson exchange approximation generating 
a Regge spectrum of relativistic bound states ^Ifl] with mass-squared 

Mli^Ml + {2n + l)\n{2 + ^/ifK^, (5) 

which is linear both in radial number n and angular momentum I. Thus the model given by Eqs.Q and Q) can 
be seen as a soluble prototype of a confining theory, here with confined fundamental fields (j) and $ and a Regge 



spectrum of relativistic bound states representing the physical particle spectrum. It is remarkable that the effective 
action for the composite fields describing the bound states can also be derived analytically within the bosonisation 
approach, thus enabling further calculations of physical quantities such as form-factors, decay widths, and scattering 
cross-sections. The more realistic choice 

1 _ p-pVa' 
Sip') ^ D{p') -5 , (6) 

where nonlocality appears as an (exponentially small at short distances) correction to the local massless propagators, 
allows a variational approximate solution of the corresponding Bethe-Salpeter equation and displays an asymptotically 
linear Regge spectrum [l^ . 

In this paper we consider the total annihilation cross section in the model of Eqs.Q and Q using the method 
suggested in |2l|, where a similar problem has been studied in application to the 't Hooft model. The total cross 
section for annihilation of a lepton-anti-lepton pair into hadrons via a scalar "photon" was evaluated in 21] by means 
of continuation of Breit-Wigner formulae in the complex plane using the fact that the physical spectrum of the 't Hooft 
model consists of a Regge spectrum of hadronic resonances with the finite widths. In our case the bound states form a 
physical Regge spectrum of resonances and the effective "hadron" action is known explicitly. The scattering problem 
can be formulated in a standard way for these hadrons and a unitary 5'-matrix can be defined. The total cross section 
is calculated by means of the optical theorem, namely via the imaginary part of the correlator of two scalar currents 

U{q) ^ J d^xe^^^O\j{x)jmO), j{x) = <i>HxMx). (7) 

As in the t'Hooft model the crucial point is to incorporate into the calculation of the RHS of Eq.© the finite widths 
of the resonances of the physical spectrum. As explained below, the imaginary part of the correlator is calculated in 
the approximation based on an extended Breit-Wigner formula. This approximation is valid if resonances are narrow 
enough, which is satisfied in the model under consideration for the seven lowest resonances corresponding to the 
energy interval s < 30 -;- 40 A^. 

It should be stressed that our task is not verification of the asymptotic inequality Eq. JQ) . At asymptotically high 
energy the number of resonances involved in the calculation of the total cross-section grows rapidly, their width 
becoming so large that they strongly overlap. In this regime an app roximation based on the Breit-Wigner formula is 
not valid and methods in the spirit of the statistical bootstrap |22l | should be more appropriate. As already stated, 
we intend rather to give a particular example demonstrating a computation of physical quantities at relatively high 
energies within a nonlocal model with confinement of fundamental fields and a Regge spectrum of bound states. 

The final result for our estimate of the spectral density is given in the upper plot of Fig[51 We conclude that in 
a nonlocal theory with confined fundamental fields and exponentially growing amplitudes, taking into account the 
physical bound states can suppress unphysical growth of observable quantities at higher energies. The finite width of 
these resonances is of particular importance. Exponential growth of the widths is a particularly distinctive feature of 
the analytic confinement scenario. It leads to the convergence of the Breit-Wigner type series for the spectral function. 
In a broader context, as indicated in a recent paper J25j, exponentially growing widths of hadronic resonances can be 
important for understanding the phenomenology related to the quark-gluon plasma: it eliminates the divergence of 
the thermodynamic functions above the Hagedorn temperature and is able to successfully describe lattice QCD data 
for the energy density. We discuss this issue briefly in the last section. 

II. THE EFFECTIVE ACTION FOR COMPOSITE FIELDS. 

Our starting point is the Euclidean functional integral 

Z[I]=M f V4>V^V^^ exp I / d^x [C{x) - ^\x)<^{x)I{x)\ j , (8) 

where the Lagrangian C{x) is defined by Eqs.Q and Q), and the external current I{x) 

I{x)^e' [ d^yG{x - y)I{y)l{y) 



represents interaction with the scalar lepton current l{y)l{y) by virtue of scalar particle (analogous to photon) exchange 
reflected in the propagator G{z). 
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FIG. 1: Graphical representation of various terms in the effective action. 

The propagators S'(p^) and D{p^) given in Eq.l^J are entire functions such that no particles can be associated 
with the fields $ and (p. These fields represent fluctuations localised in space and time. A typical space-time size of 
fluctuations is set by a confinement scale A. The physical particle spectrum of the system can be recovered by means 
of the "bosonisation" procedure applied to the functional integral Eq. ((SJ . For a detailed description of this procedure 
for this case we refer the reader to |l9f . Essentially, bosonisation requires solving the Bethe-Salpeter equation in the 
one-boson exchange approximation, which can be done analytically due to the simple Gaussian form of propagators. 
The main steps are as follows. Integrating out the field (j) one arrives at a quartic interaction of the scalar fields <&: 

L2m = Y f d^xid''x2^\xiMxi)D{xi - a;2)$t(a;2)$(a;2). (9) 

Introducing a complete orthonormal set of functions Uq, corresponding to the radial quantum number, total momen- 
tum and magnetic numbers Q = {n,l,fj,}, the four point interaction can be rewritten as an infinite sum of products 
of a non-local currents 

L2 = ^5]/d'xJQ(a;)jQ(a:) (10) 

with 

Jq{x) = $t(x)T/s(a)$(a;) 



Vaid) = / d^yVD{y)UQ{y)e^^ (11) 



One introduces auxiliary fields $q with a subsequent Gaussian integration over the field $. As a result the original 
functional integral IjHJl turns out to be identically rewritten in terms of composite fields ^g representing collective 
excitations in the system 

where the dimensionless coupling constant a — (7^/(47rA)^, and 

Wi[g^] = -Tr[ln(l - g^aVaS - IS) + ^^ qVqS^ q,Vq> S] (13) 

describes interactions between composite fields from which the quadratic term in the expansion of the logarithm is 
subtracted and added to the original quadratic part, 

aI]QQ'(x-2/) = ^j^VQ{d^)S{x-y)VQ'{dy)S{y-x) 



The nontrivial requirement is that the basis Uq be chosen such that the self-energy T,qqi (p) and hence the quadratic 
part of the effective action is diagonal in the quantum numbers Q, 



^QQ'ip) ^ Eq{-p^)5qq. 



Bound state masses are then real solutions to 



l = aEQ{Ml). 



(14) 



(15) 



Diagonalisation of the self-energy is equivalent to the solution of the Bethe-Salpeter equation in the one-boson exchange 
approximation. In the momentum representation propagators and vertices involved in the effective action have the 
form [IJj 
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FIG. 2: The two-particle decay width r„(M^) as a function of radial quantum number. Deviation from strictly exponential 
behaviour is due to the polynomial Pn{s) in r„(s) oc P„(s) exp(const s). 



S{p) = A-2e-pVA^ 
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{-ir+^kC^iTl^_^^{K)L^+\aK')e 



'bK^ 



(16) 



with K = {p + q)/2A, a = 2^3, b = 2/(1 + V3) = 4/(2 + a) and 



(^n 



n\{l + l) 



(l + V3)2"+'+2 ^ {n + l + l)\' 
Here the angular part of the vertex is given by T' , the irreducible tensors of the Euclidean rotation group 0(4), 
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with C; being Gegenbauer polynomials. The radial part corresponds to the Laguerre polynomials L^^^, with 



duu^+h'''Ll+\u)Ll+\u) = S. 



Calculation of the self-energy gives 



Enl{-P^) = 



e 2A-i 



(2 + V3)2"+'+2 ' 



and the square of the bound state masses read 



M?,i = 2A2 



In ^^ + ^^ + {2n + I) ln(2 + VS) 



(18) 



which manifests a linear Regge spectrum. Finally the fields should be rescaled 

A2 di?Q(-p2) 



h~^ - 
S " (4^)2 



dp2 



p2 = ^M2 



(19) 



in order to ensure the correct residue of the propagator at the mass pole. The coupling constant /ig is defined as a 
dimensionless quantity. Thus the final form of the functional integral for composite fields is 



^ = li I T^^^Q^^pI^^^ J d^P'^Qi-p)i^-^EQ{-p^))'i'Q{p) + Wi[hQ^ 



(20) 



The original coupling constant g enters only the quadratic part of the effective action, while the remaining terms 
contain the effective coupling constant Hq. 

The formal functional integral H2U|) can be used to define a unitary nonlocal theory for the interacting fields ^g by 
defining the appropriate Gaussian measure for its computation. Namely, one rewrites Eq{—p'^) as 



1 - aEQi^p^) - {Ml+p^)aE'Q{M^Q) - [aSg(V) - 1 + {M^q + p'')aE'Q{Ml)] 
1 - aEQ{Ml) 



(21) 



0. 



This decomposition exists since Eq{~p'^) is an entire analytical function. If one uses now the first term of the RHS of 
Eq. (|21|1 to define the standard free field Gaussian measure for computation of the functional integral then the theory of 
interacting composite fields ^g so defined is unitary since all nonlocal form factors in the action are entire analytical 
functions, and the only singularities which appear in the S'-matrix are related to the physical particles associated with 
fields "fg. In particular, this theory satisfies the conditions required for the derivation of the high energy bound Eq.fQ) 
with K — 1. Such a "free field" Gaussian measure is appropriate for computing the processes in which all composite 
particles are almost on-shell. In this case [a-Bg(— p^) — 1 + {Mq + p'^)aE'Q{MQ)] = 0{{Mq -fp^)^) is small and can 
be considered as a perturbation. 
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FIG. 3: Graphical representation of the polarisation function. 



III. TOTAL CROSS SECTION OF // ANNIHILATION TO "HADRONS". 



According to the optical theorem the total cross section for the inclusive process II -^ "hadrons" 



can be expressed through the imaginary part of the amplitude for elastic forward scattering, which is proportional to 
the spectral density p{s), namely the imaginary part of the correlator of two scalar currents 



p{s)=-s d'xe^'>^0\j{x)jmO)\ 



In terms of the effective action in Ea. H12l) the spectral density is given by the imaginary part of diagrams with one, 
two and generally n intermediate "hadrons" as shown in Fig|3Ji-c. Since the propagator of the constituent field <& 
(thin solid line) and the vertices Vq (see Ea. ()16|) ) are entire analytical functions no contribution to the imaginary part 
of diagrams in Fig|31can come from the constituent field loops. In particular, the diagram in FiglSJl has no imaginary 
part and does not contribute to the spectral density. 




FIG. 4: Spectral density in the Breit-Wigner approximation as defined by Eg. 1)22^ 

In the Breit-Wigner approximation of the resonance propagator, the contribution of diagram in FiglSt^ to the 
spectral density takes the form 



.^ ^1- fenOg^o(Mno)[£ + r„(M^o)] 



(22) 
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where r„(M^Q) is the width of the resonance, which is equal to zero for the lowest state. It should be noted that the 
expansion analogous to that in l/N^ in the 't Hooft model or QCD is represented in the model under consideration 
via an expansion in the number of loops which include internal resonance propagators. Further calculations will be 
based on this decomposition over "hadronic" loops. At this point we take this expansion as a way of classifying the 
various contributions, and will discuss its validity in more detail shortly. As for the leading order contributions in 
1/A^c in the 't Hooft model, the bound states in the nonlocal model are stable in the zeroeth order of this expansion 
over "hadronic" loops. In this approximation, namely with decay widths neglected in Ea. l|22|l . the spectral density is 
simply an infinite equidistant sequence of delta functions. The lowest order contribution to the widths of resonances 
in this loop expansion are the two-particle decays. Numerically the widths for few lowest resonances are shown in 
FigEl One can see that the width is exponentially growing with n, which ensures convergence of the series for the 
spectral density in Ea. l|^ for any s since 

l2 2 
"-noSnO °^ '* 

for n 3> 1. The exponential growth of the decay width is precisely due to nonlocality. This should be compared with 
the behaviour r„ oc \/n and gn — *■ const in the 't Hooft model or string based picture of confinement. We comment in 
the discussion section on whether available experimental data can distinguish between these two qualitatively different 
dependencies. 

Summing the convergent series in Ea. (|22|l gives the spectral density shown in Fig^ Qualitatively it has a reasonable 
form and decreases at large s. The advantage of the Breit-Wigner approximation for the resonance propagator is that 
it explicitly ensures the correct analytical properties: only physical singularities in the complex s-plane are present in 
the Breit-Wigner propagator. 

However, there are three important factors at large s which are missed in the spectral density calculated through the 
Breit-Wigner propagator with inclusion of two-particle decay widths. First of all decays to three and more particles 
become important at large s. Accounting for them would broaden the resonances and lower the spectral density. 
The second contribution comes from the diagrams with two or more intermediate bound states (see Fig|3|D,c) which 
would increase the spectral density at large s. These two factors indicate that Ea. (|22|l cannot give a truly reliable 
result for asymptotically large s. At most it provides a hint at the tendency of p(s) to decrease at larger s. The third 
factor is that the Breit-Wigner approximation is valid, strictly speaking, only in the vicinity of a resonance pole and 
if the resonance is narrow enough. This means for s located away from the resonance position the contribution of the 
n-ih resonance to the series Eq.(j221) might not be reflected correctly. In terms of the effective action discussed in the 
previous section the Breit-Wigner approximation for the propagator corresponds to the "free field" Gaussian measure 
in the functional integral with the addition of the on-shell widths of resonances. It is clear that corrections coming 
from the second term in the RHS of Eq.|(2U arc not small for p^ being off-shell. Certain resummations are necessary. 

The effective action at our disposal offers an extension of the Breit-Wigner representation, which can account for 
such resummations just mentioned via inclusion of the energy dependence of the self-energy £'„o(s), the transition 
amplitude (?no(s) and the imaginary part - decay width r„o(s). This extension of the spectral density takes the form 

oo 

P{s) = ^Pn{s), 

n=0 

|5no(5)P[£-t-f„o(s)] 



Pn{s) = lim 



^nO = 7 — ^^^"O (23) 

where the summation spans only radial excitations with I = because of conservation of orbital momentum which is 
a good quantum number in this problem. The transition amplitude gnois) is given by the diagram in Fig^a. Simple 
calculation gives 

gs/2A2 

9no{s) = i-irV^^TlE^ois), EMs) = ^^ „^, . (24) 

The quantity r„o(s) is the total width of two-particle decays of a resonance with orbital momentum I = 0, radial 
number n and squared-energy s. Details of the calculation of r„o(s) are given in the appendix. 
Numerical results for r„o(s) are given in FigEl which shows its exponential growth with s, and 

r„o(s) > ff«o(s) oc Eno{s) 
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FIG. 5: Spectral density p(s) as a function of the energy variable s (upper plot) and the contributions pn{s) of few lowest radial 
numbers to the total spectral density (lower plot) 



for large s. This inequality results in the behaviour of the separate terms Pn(s) shown in the lower plot in Fig|Sl 
Summing up these terms gives the spectral density shown in the upper plot in Fig|31 Qualitatively this spectral 
density reproduces the character of the simple Breit-Wigner approximation, FigQ] The oscillations seen in Pn{s) in 
Figl^l additional to the resonance peaks are reflections of two features. The first is the corresponding oscillations in 
r„(s) in Figl^ which are due to the opening of new decay channels at certain values of s. The second, and more 
important, cause is the functional structure of the amplitude An{s) (see appendix): essentially it is a product of a 
polynomial and exponent in s. Again the spectral density shows several sharp peaks corresponding to the lowest 
narrow resonances, then it has a maximum and a tendency to decrease at higher energy. It should be stressed that 
/9(s) is limited by the region of not very large values of s because of the same as for pbwCs) two reasons (decays into 
more than two final states and diagrams with two or more intermediate resonances). However unlike the plain Breit- 
Wigner approximation the extended formulae for the propagator should be more appropriate for the values of s away 
from resonance position although only in the vicinity of the real axis in the complex s-plane. The reason preventing 
it from being a good approximation in the whole complex plane is hidden in the property of entire functions which 
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FIG. 6: Two-particle decay widths rn(s) as a function of energy variable s for the few lowest resonances. 

have infinitely many zeroes in the complex plane. In addition to the desirable physical pole corresponding to the 
resonance, the extended propagator has artificial unphysical poles for complex values of s. As soon as these poles are 
sufhciently far from the real axis and do not noticeably affect the form of p{s) the approximation Eq. I|23|) can be 
considered appropriate. In the absence of the term with r„ in Eq. if^ and using Eq. ((TH|l . the unphysical singularities 
would appear at s = M^^^ + Aiirkh? for any nonzero integer k. The addition of r„(s) shifts all poles to the complex 
plane, and a cursory numerical check reveals that the artificial singularities produced by the resummation are located 
sufficiently far from the real axis. 



IV. DISCUSSION 



In summary, we have presented the spectral density in the Breit-Wigner approximation for a nonlocal model with 
confined constituent fields but a Regge spectrum of bound states. The two-particle decay widths have been computed 
and incorporated into the evaluation of the spectral density. The main observations are: firstly, due to rapid growth of 
resonance widths with growing radial number the formal sum over resonances (which represents the spectral density) is 
convergent; secondly, apart from the expected delta- function like peaks corresponding to the lowest narrow resonances 
the spectral density in this approximation decays at larger values of the energy variable s as shown in Fig^ 

The Breit-Wigner approximation roughly reproduces the behaviour of the resonance propagator as a function of 
s away from the resonance position. In order to improve this a resummation of the series over the effective "bound 
state-constituent field" coupling constants ft-g are necessary, what we have called an "extended" Breit-Wigner ap- 
proximation. The spectral density computed by means of such a resummation shows qualitatively the same behaviour 
as the Breit-Wigner spectral function - there are several peaks corresponding to the lowest narrow resonances and a 
tendency to decay at higher energies. 

As already discussed, the resummation introduces artificial unphysical singularities. However, in the particular 
case under consideration these unphysical singularities are located sufficiently far from the real s-axis. Insofar as we 
consider the resonance propagator in the vicinity of the real axis, the approximation based on resummation of higher 



11 

orders in Hq is quite reliable. It should be noted that, in the context of nonlocal quantum field theory, the problem of 
unphysical singularities appearing in the complex s-plane is akin to the Landau pole problem in local quantum field 
theory: a partial resummation of the perturbation series is generally hard to achieve without introducing unphysical 
singularities, which applies both to local and nonlocal models. There are two aspects to this issue. The practical 
way to manage is to look for summation prescriptions which are appropriate in the restricted region of the complex 
momentum plane relevant to a given problem; this is what we have pursued in the present work. A more fundamental 
approach to the problem would require formulation of the general principles and methods of summation which would 
avoid the appearance of unphysical singularities. 

As mentioned, the asymptotic dependence of two particle decay widths on the bound state quantum number in the 
nonlocal model r„ ex expjconst n} is drastically different from the dependence in the t'Hooft model r„ ex ^/n, and 
semi-classical string based estimations of 23j where T{M) oc M. The computation in the 't Hooft model is similar to 
ours in that decays of excitations to lower states at the same trajectory are considered, while the computation of |23l | 
applies to all possible decays. It is interesting to take a look at experimental data on decay widths in the hope to 
distinguish between these qualitatively different situations. Unfortunately the available data seem to be inconclusive 
in this respect. The problem is that for the purpose of extracting asymptotic dependence of the width on radial or 
orbital number one should consider only decay modes of a given excited state to states lower on the same Regge 
trajectory. In the real multiflavour world, decays of excited states are dominated by the modes into final states which 
include ground state mesons from Regge trajectories of other flavour octet states. A typical example is the family of 
orbital excitations of K* , 

{{K* {892), k;{U30),k; (1780), Kl{2045)) 

with the total widths (using [23| and data available through the 2003 partial update) 

Tk' = 50 ± .9MeV, Ta- = 98.5 ± 2.9MeV, Tk^ = 159 ± .21MeV, T^* = 198 ± 30MeV, 

looking as well fitting linear in I behaviour of the width (though T ex e*^ expected in the domain model [14 . Il5j | would 
also fit but be implausible with the limited range of data). However the problem is that decay modes which contribute 
to the widths of these resonances all include K, tt, p, i] mesons which obviously do not belong to the trajectory of 
K*. So, this approximately linear dependence in I does not contain the information which we are interested in. 
It can be strongly infiuenced by chiral dynamics, and its description requires a model with both chiral symmetry 
and confinement implemented simultaneously (for instance the model of 1^ has a chance to provide such a suitable 
framework for calculating partial widths known from experiment and enabling thus a detailed analysis). The clearest 
case would be the Regge family of radial excitations of a pseudoscalar meson like tt, 7r(1300),7r(1800), however the 
widths and especially the partial widths for different decay modes are not known for these mesons with the required 
accuracy. These two examples are quite typical, and we conclude that it is not yet possible to distinguish between 
the two asymptotic behaviours on the basis of experimental data. It should be stressed that the nonlocal model, 
despite the exponential asymptotics for the width, indicates very small partial widths of the few lowest resonances 
for decays into lower states on the same trajectory due to the polynomial pre-exponent, and can considerably deviate 
from strictly exponential form as can be seen from Fig. These partial modes would anyway be screened by other 
modes in the multifiavour real world. A further signature of exponential behaviour for larger (n, I) states would 
be that observable trajectories cannot be "long", since resonances with radial or orbital numbers bigger than some 
critical ric {ric = 6 in Fig. |^ or Ic are too wide to be observable. Such a sharp cut-off of the trajectories would be 
typical for exponential dependence, but not for linear or square root behaviour. The property that the known Regge 
trajectories, for example, for unflavoured mesons cut-off at about 2 GeV may be an indication of this feature. 

The crucial role of the finite widths in modulating the behaviour of the cross-section for growing energy can be 
understood in terms of the energy being dissipated into the creation of unstable particles, which decay into more 
stable particles and so on. In this respect, the mechanism is close in spirit to the Hagedorn mechanism for a 
maximal temperature in hadron-hadron collisions |22l | where one has no recourse to the optical theorem and thus the 
machinery of the statistical bootstrap is necessary a priori. In our case, this framework would be useful once the 
narrow width approximation has broken down. In this context we mention that exponentially growing decay widths 
can impact on the deconfinement phase transition in particular the persistence of hadronic properties beyond the 
critical temperature. For example in |25|, an ansatz involving exp{rn/TH) for the width contribution to the spectral 
function of a hadron gas at temperatures greater than the Hagedorn temperature, Th, reproduces the characteristic 
slow approach to the Stefan-Boltzmann ideal quark-gluon gas result seen in lattice calculations 26). Interestingly, the 
Hagedorn temperature scale, traditionally argued to be determined by the mass of the lightest hadron, is commensurate 
with the confinement scale which figures prominently in analytical confinement and the exponentially growing decay 
widths discussed in this paper: A ^ Th ^ 200MeV. The possibility that the confinement scale is more significant in 
fixing the Hagedorn temperature, especially in the chiral limit, has been argued recently [27] . Further study of the 
consequences of analytic confinement to high temperature behaviour is appropriate. 
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Finally, we have not addressed the relation between our study to the idea and formalism of "quark-hadron duality" . 
With the purely Gaussian propagators treated in this article, such an investigation would be pointless. However this 
problem becomes well- formulated if one considers the more realistic propagators of the form Eq. ^ . This propagator 
has the same short distance behaviour as the local one and nonlocal corrections are exponentially small in the deep 
Euclidean region. With such propagators one can address the issue of q uark-hadron duality within a model with 
dynamical confinement of constituent fields, in a spirit close to that of 



n q ui 



APPENDIX A: TWO-PARTICLE DECAY WIDTH. 

The width for two-particle decay for a state with the time-like four momentum squared p^ — s decaying to particles 
with masses rrii and m,2 is given by the standard equation 



T(^)[s,m,,m^) = ^^^j n '^/'\p-Pi-P2)\A{p,pr,p2)\''. (Al) 



The modulus of the amplitude is independent of the momenta. Thus for the width for the decay of the state with 
angular momentum I = and radial number n into the final states with masses rrii and r7i2 and quantum numbers 
(ni,Zi) and (712,^2) respectively, we get 

p(2) _ 1 R2{s,mi,m2) ^ , 2 2x 

'- n,nih,n2l2 ~ (n-w']'^ 2s -^ n,niii ,n2i2 W; "^1 7 "^2^ 

TnMih,n2l2 = |^n,riih,«2i2(PiPliP2)| , 

where i?2(s, 7711,7712) is the two-particle phase space 

-R2(s, 7771,7712) = —A2(s, nil, 7712) 



The total width is then 



A(s, 7/11,7712) — {s — {mi + 77i2)^)(s — (?7ii — 7712)^). 

ni/i, 712/1 

The amplitude An-^mh. 11212 which should be computed is given by the diagram Fig^. For arbitrary (li^h) and 
polarisations, the corresponding loop integral takes the form (all momenta below represent dimensionless ratios, e.g. 
P = P/A) 

^;^;".^;C-i'....,, = ^KoK,lK2l J ■^V""{q)S{q+p/2)V;:i%^^ {q+p/2 - h/2)S{q + p/2 - h) 

xV:[^X^q-h/2)S{q-p/2). 

Since the orbital momentum of the incoming state is zero, the final states are required to be on-shell, and the angular 
momentum is a conserved quantum number in this problem, we need to extract from this amplitude the contribution 
of final states with total angular momentum equal to zero. Thus we need the term with the angular momentum I — 
in the decomposition of the amplitude ^^^^^f^'^y'l y^ over the irreducible tensors T' with I = 0, . . . ,Zi + ^2- The 
required contribution is proportional to the trace of ^ui"^uV"t'i^ vi f^^' ^2 = ^1 

A-l 

A , _ \ ^ An\nxl\,n2l\ 

■^n;ni,n2;h — ,, ^ ^^^l ■■■A'!i .Pi •••A'ii ' 
ti . 

Ml-'-Mii 

where A is defined as a dimensionless quantity. Taking into account Eqs.lQSJ and H17(l we arrive at 

An;nun2;l = TJ^j KoKilhn2lGnoCmlCn2l / T^^Q [{^ + {P ^ ki) /2){q ~ ki /2)] 

X Ll[aq^]L\tM'i +iP- fci)/2)^]i'+M«(9 - h/2)']Ci'\{q + {p - fci)/2)(q - fci/2)] 

X exp {-bq^ - b{q +{p- fci)/2)2 - biq - ki/2f - {q+p/2f - {q + p/2 - k^f - (q - p/2f] , (A2) 
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where C; are Gegenbauer polynomials. The Gaussian integral over the loop momentum q can be performed analyt- 
ically. The most convenient way to do this is to use a generating function for Laguerre polynomials 



E 

n=0 



L^{x)z -- ^ _ z)'3+i' 



Both sides of Ea. (|A2|l are multiplied by z"z"^Z2^ and summed over n, ni, n2 and the result of integration is expanded 
in a power series in each of z, zi and Z2- The coefficients of expansion gives the amplitudes An:n-i^,n2\i- The Gegenbauer 
polynomials are substituted in explicit form. This procedure can be easily implemented in Mathematica or Maple. 
Finally the amplitude An-ni,n2;i 3.S a function of s is a product of a polynomial in s and an exponential of s. The 
degree of the polynomial depends on n, ni, 712 and I. 

The total two-particle decay width of the resonance with four momentum squared s, zero orbital momentum and 
radial number n entering Ea. (|23|l is then 

A^ 

/,ni ,n2 

where all states with nonzero phase space are summed. 
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